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Abstract. We report on the detection of a rapidly variable narrow Fe Ka line in Mkn 841. The source has been 
observed two times by XMM-Newton and simultaneously with BeppoSAX. The two observations, of about lOks 
long each, were separated by ~ 15 hours. The line flux reaches a maximum during the first observation and is 
significantly reduced in the second one. The continuum shape and flux, instead, keep roughly constant between 
the two pointings. Such rapid variability of a narrow (unresolved by the XMM-pn instrument) line has never been 
reported in the past. These results are not easily explained in the standard cold reflection model where the narrow 
line component is supposed to be produced far from the primary X-ray source (e.g. from the torus) and is thus 
not expected to vary rapidly. Different interpretations are discussed. 
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1. Introduction 

The fluorescent Ka iron line is an important feature in the 
high energy spectra of Seyfert galaxies. It is believed to be 
the result of fluorescent emission due to photoionization of 
iron atoms in optically thick matter illuminated by a com- 
pact source of hard X-ray photons. The study of the iron 
line properties (ionization state, profile, variability) indeed 
provides very useful information on the geometry and the 
physical properties of the emitting material (Fabian et al. 
200C and references therein). 

Rapid variability (on less than a day time scale) of the 
iron line has been reported in the past for a few objects 
(Yaqoob et al. 1996; Vaughan & Edelson 2001: Wang et al. 
2001). In all these cases, the line was relatively broad. This 
broadening is believed to be a direct probe of Doppler and 
gravitational effects that occur in the region very close to 
the putative supermassive black hole (hereafter BH) . The 
line variability would thus be associated with the strong 
perturbations of the emitting medium in the vicinity of the 
BH. Moreover, a narrow iron line appears to be often, if 
not always, present in Seyfert 1 galaxies (Pounds & Reeves 
2002 and references therein). It is thought to be produced 



by an obscuring and remote torus postulated to exist in 



AGN unification schemes (Antonucci [1993 ; Ghisellini et 
al. 1994). This line is thus not expected to vary on so 
short time scale. 

Mkn 841 is a bright Seyfert 1 galaxy (z=0.0365), one of 
the rare Seyfert 1 detected by OSSE at more than 3 a 



(Johnson |1997 ). It is known for its la rge sp ectral variabil 



ity (George et al. |l993t Nandra et al. 1995 ), its strong soft 



excess (this was the first obje ct where a soft excess was 



observed, Arnaud et al. 1985) 



and its variable iron line 
The 



1993) 



(at least on a year time scale, George et al 
latter was observed in some cases with a relatively large 
eq uivale nt width (hereafter EW) of about 400 eV (Day et 
al " 



1990 



Bianchi et al. 2001 ) significantly above the value 
predicted by st andard cold reflection model (e.g. George 
& Fabian 



1991|) 



In this paper we present the results from a XMM-Newton 
observation of Mkn 841. This source was also observed si- 
multaneously with the BeppoSAX satellite. We report the 
presence of a rapidly variable narrow (i.e. unresolved by 
the XMM-pn instrument) iron line changing by a factor 
of a few in flux on a half day time scale. 
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MECS 2-10 keV light curve 
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Fig. 1. MECS light curves (bins of 2500 s). The periods 
of the two XMM observations are indicated on the figure. 



two simultaneous XMM observations are also reported in 
this figure. The EPIC-pn light curves show smooth flux 
variations of a few percent on time scale of 1000 s or less. 
The integrated XMM-Newton flux of Mkn 841 over the 2- 
10 keV range was roughly the same for the two EPIC-pn 
observations at about 1.4xl0~ n erg. s _1 . cm -2 . It slightly 
differs by ~10% (smaller) than the mean 2-10 keV flux 
measured in the MECS as it is generally observed between 
these two instruments (Molendi private communication). 
The (2-5 keV)/(5-10 keV) hardness ratio of the MECS 
being consistent with no spectral variability during the 
observation, and in order to maximize the statistics, we 
integrated the MECS and PDS spectrum over the entire 
observation despite only a part of it was simultaneous with 
the XMM one. In the following, all errors refer to 90% con- 
fidence level for 1 interesting parameter (A% 2 =2.7). 



2. Observation and data analysis 

Due to operational contingency, the 30ks observation of 
Mkn 841 with XMM (Jansen et al. |200l| and references 
therein) was split into two parts, planned the 13th (noted 
OBSl) and the 14th (noted OBS2) of February 2001 with 
~ 12 and 15 ks exposure time respectively. The two ob- 
servations were separated by about 15 hours. 
In this paper, we will only deal with the EPIC-pn data (we 
have checked that the pn and MOS spectra are consistent 
with each other in the 2-10 keV band). A complete anal- 
ysis, including all the XMM instruments, is deferred to a 
future work. The EPIC-pn camera was operated in Small 
Window mode (hereafter SW) , with thin aluminium filters 
to reject visible light. The EPIC-pn event file was repro- 
cessed from the ODF data files using the epchain pipeline 
task of the XMM Science Analysis System (SASv5.2) and 
using the most updated version of the public calibration 
files. We note that SASv5.2 only contained a prelimi- 
nary model for the pn-chips charge transfer efficiency. This 
should not affect our main results presented here, based 
on the difference between two spectra, but may change 
slightly the spectral shape. The source spectra and light 
curves were built from photons detected within a 40 arcsec 
extraction window centered on the source. X-ray events 
corresponding to pattern < 4 were selected. The back- 
ground was estimated from an offset position and was 
found to be extremely low during the whole observation. 
The EPIC-pn count rate was of about 18 cts.s -1 for both 
observations, i.e. well below the pile-up threshold of ~ 300 



3. The line variability 

The XMM-Newton and BeppoSAX instruments need to 
be carefully normalized one with each other to correct 
from calibration uncertainties. The MECS data were used 
to perform the cross calibration as follows. The multi- 
plicative factor between the EPIC-pn and the PDS was 
obtained by normalizing the EPIC-pn flux to the MECS 
spectrum and by applying a relative normalization factor 
of 0.86 (free to vary by ±5%) corresponding to the cross 
calibration of the MECS and PDS instruments (Fiore et 
al. 1999). We then obtain a normalization factor of 1.13 



and 1.08 for the PDS spectrum relative to the EPIC-pn 
one for OBSl and OBS2 respectively. From now, we will 
only use the EPIC-pn and PDS instruments after having 
checked that the EPIC-pn and MECS spectral shape were 
consistent with each other in the 2-10 keV band. 
We have reported in Fig. || the ratio between the EPIC- 
pn spectra of OBSl and OBS2 and the best fit pexrav 
model of xspec (Magdziarz & Zdziarski 1995 ). This best 
fit was obtained fitting the EPIC-pn and PDS data simul- 
taneously. A line is clearly visible near 6.2 keV (~ 6.4 keV 
in source rest frame) in OBSl. Indeed the addition of a 
Gaussian in the fit of OBS 1 is significant at more than 
99% (following the F-test) with Ax 2 =20 for 3 additional 
parameters. On the contrary, the addition of a line was 
not significant in OBS2 (Ax 2 = l). The best fit parameters 
of the PEXRAV+GAUSSIAN fits are reported in Table |l|. 
The energy of the line detected in OBSl is consistent with 
a neutral Ka line £^ me =6.41^Q'gg keV (source frame), it is 



alent width (EW) o f 120 +^ eV (fitting with a diskline 



cts.s -1 in SW mode. For the fitting procedure the XMM unresolved with a width cr lme <170 eV and has an equiv- 
spectra were binned to have at least 50 counts per bin. 
Mkn 841 was also observed between the 11th and 14th 
of February 2001 by the BeppoSAX satellite (Boella et 
al. 1997 and references therein). The source was pointed 
for a total net time exposure of ~90, 40 and 20 ks for the 
MECS, PDS and LECS instruments respectively. Since we 
focus our analysis, in this paper, to data above 2 keV, we 
will not deal with the LECS data any longer. The MECS 
light curve is reported in Fig. [l] with a binning of 2500 s. 
Variations of about 20% are observed. The periods of the 



model (Fabian et al. 1989), the data still constrain the line 
to have a narrow profile, the inner disk radius being larger 
than —200 Schwarzschild radii). We obtain an upper limit 
of -70 eV for the line EW in OBS2 and the best fit line 
flux is a factor ~3 weaker than in OBSl. 
We then checked for line flux variations on smaller time 
scales. For that purpose, we studied the light curves of 
the flux ratios (6-6.6 keV)/(3-5 keV) and (6-6.6 keV)/(7- 
10 keV). These ratios are expected to be sensitive to the 



P.O. Petrucci et al.: A rapidly variable narrow X-ray iron line in Mkn 841 



3 



' * " " I 

i 

■ah, i 




XMM Observation 1 


i 





"> I , , , , i , , , , i , , , , I 

°5 6 7 8 



Energy (keV) 





i i i i i | i i i i i i i i i | i i i i 
XMM Observation 2 







° 5 5.5 6 6.5 7 7.5 8 

Energy (keV) 



Fig. 2. Ratios data/model (without line) for OBS1 and 
OBS2. The PEXRAV fit was done using the EPIC-pn and 
PDS data simultaneously. The vertical dashed line indi- 
cates the position of the 6.4 keV neutral iron line in the 
source rest frame. 
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Table 1. Best fit parameters for OBS1 and OBS2 when 
fitting with the PEXRAV+GAUSSIAN model above 2 keV. 
The line energy is in the source rest frame. 



line flux variations with respect to the continuum shape 
in each side of the line. The line interval is smaller than 
the others to be sensible enough to the line flux variability. 
The light curves show a significant increase of both flux ra- 
tios during OBS1 where the line flux reaches a maximum 
of ~5xl(T 5 ph.cm -2 .s _:t . Fits of the different EPIC-pn 
spectra along the light curves (using a temporal bining of 
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Fig. 3. Contour plot (68%, 90% and 99% confidence level) 
of the line flux vs. the line energy during the line flux 
maximum of OBS1 (solid line) and the total OBS2 (dashed 
line). Cf. text for details. 

1800 s) show that the best fit line flux varies by a factor 
~ 6 during the total observation (the lowest flux being in 
OBS2). The line width is always consistent with zero i.e. 
the line is always unresolved by the instrument. 
We have plotted on Fig. || the contour plots of the line 
flux versus the line energy (with contours at 68, 90 and 
99%) obtained during the sub-interval of OBS1 where the 
line flux reaches its maximum (solid contours) and for the 
total OBS2 (dashed contours). These contour plots have 
been obtained using the EPIC-pn data alone fitted with 
PEXRAV + GAUSSIAN. The line width was kept fixed to 
the value of 0.1 keV and we also fixed the reflection nor- 
malization parameter R to its best fit value. From these 
contours, the hypothesis of the constancy of the line can be 
rejected at more than 99.9% (it is rejected at 99.5% if we 
use the total OBS1). Unfortunately, the BeppoSAX/PDS 
data does not provide enough counts to test whether the 
reflection also varied significantly during the pointing. 

4. Discussion 

Two observational facts of the present data are not easy 
to account for. Firstly the narrow line varies rapidly by 
a factor of a few while the continuum changes by only 
~ 10-20%. Secondly the reflection seems relatively large, 
R > 1, while the line EW is relatively small < 170 eV 
(it is worth noting that fitting the BeppoSAX data above 
2 keV gives also a large R=2.bt_\'\ and a small EW < 
100 eV). These points are in contradiction with the stan- 
dard cold reflection model. In this model, the line width 
is due to the motion of the emitting matter in the BH 
gravitational potential. In this case, a fluorescent neutral 
iron line with a width of about 150 eV (~ 14000 km.s -1 ), 
like in our case, is expected to be emitted at a distance 
r ~ 10 17 M 8 /q cm where the BH mass Af bn = 10 8 M 8 so- 
lar masses and q depends on the distribution of the orbital 
shape and inclination of the matter spiraling around the 
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BH (q is expected to be larger than 1, Krolik 2001). 
On the other hand, a variability time scale of ~ 10 hours 
constrains the size d of the line emitting region to about 10 
light hours, i.e. d ~ 10 15 cm. The solid angle sustained by 
the emitting region, as seen by the X-ray source, is then 
AO. ~ Trd 2 /r 2 ~ 3 x 10~V^V 2 - The observed values 
of R = Afl/2ir >1 then imply gM 8 " 1 >140. However, the 
fact that we do not observe any flare in the MECS light 
curve, which would be at the origin of the line variability, 
constrains r to be larger than about 1 light day, meaning 
qM% l <4Q in contradiction with the previous estimation 
and independently of the BH mass. Modifications of this 
standard view are thus clearly needed. 
It is first possible that the observed line variability is not 
controlled by the continuum itself but by an intrinsic vari- 
ability of the reflecting material due to, e.g. thermal insta- 
bilities. These instabilities may change the ionization state 
of the upper layers of the matter. The strength of the neu- 
tral iron line, expected to be produced in the deeper layers, 
could then be degraded by th e skin scattering prop erties 
(Nayakshin & Kallman |200l[ Petrucci et a l. |200l|) . On 



the other hand, Nayakshin & Kazanas ( 2002) ) also showed 
that, in response to changes in the X-ray flux, the time 
for re-adjustment of the hydrostatic balance of the illu- 
minated gas may be longer than the light-travel time be- 
tween the primary X-ray source and the reflecting mate- 
rial. The line variability observed here may then result 
from an X-ray flare produced in the inner disc region 
but occuring some time before the start of the XMM and 
even BeppoSAX observations. We note however that the 
change of ionization state of the emitting region suggested 
above would probably produce a ionized iron line which is 
not detected. 

Alternatively, the line emission may have been affected 
by the continuum outside of the observed band pass. 
Variations of the high-energy cut-off and/or of the UV- 
EUV emission can indeed alter the ionization structure of 
the disc surface without obvious chang es in the 2-10 keV 
continuum (Vaughan & Edelson |200lf) . The EPIC data 
down to ^0.5 keV do not show significant variation be- 
tween OBS1 and OBS2 but a more precise study of the 
RGS and OM data will help to better test this hypothesis. 
Then, if we admit that the line variability can exist with- 
out significant variations of the continuum in the XMM 
energy range, the distance r of the emitting line region 
estimated above could be compatible with a half day vari- 
ability time scale (or less) if the BH mass of Mkn 841 
is relatively small, < 10 6 solar masses. It is however sig- 
nificantly smaller than the ~ 10 s solar masses recently 
estimated by Laor (2001). 

Interestingly, a narrow profile is also expected, indepen- 
dently of the BH mass, if the flare at the origin of the 
line variability occurs in the inner regions but close to the 
disc and thus illuminates a small fraction of its surface 
(Nayakshin & Kazanas |2001 ; Yaqoob 2001). However, in 
this case the line is generally expected to be shifted due 
to relativistic effects, which is not observed. 
A less "standard" explanation could be the crossing of an 



obscuring thick cloud in between the X-ray source and the 
reflecting matter. This cloud needs to be close enough to 
the X-ray source to rapidly eclipse a large part of its emis- 
sion. It can explain why the continuum does not change 
between the two XMM observations since the cloud may 
not cross the line of sight. We then would expect the reflec- 
tion hump to vary, an hypothesis which is unfortunately 
not testable with the present data. 

5. Conclusion 

The neutral narrow iron line of Mkn 841 appears to be 
highly variable on a timescale of 10 hours or less. This 
cannot be easily explained in the standard cold reflection 
model framework. None of the explanations proposed here 
appear completely conclusive. A consistent interpretation 
of these results has to concile the presence of this variable 
narrow and apparently neutral line, a constant (in flux 
and shape) underlying X-ray continuum and a large (>1) 
reflection normalization. New data are clearly needed to 
confirm this observation and to understand the origin of 
the line variability. 

Acknowledgements. We thanks the referee, S. Vaughan, for his 
useful comments. 
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